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Abstract 

The local correlation between far-infrared (FIR) 
emission and radio-continuum (RC) emission for the 
Small Magellanic Cloud (SMC) is investigated over 
scales from 3 kpc to 0.01 kpc. Here, we report good 
FIR/RC correlation down to ~15 pc. The reciprocal 
slope of the FIR/RC emission correlation (RC/FIR) in 
the SMC is shown to be greatest in the most active star 
forming regions with a power law slope of ^1.14 indi- 
cating that the RC emission increases faster than the 
FIR emission. The slope of the other regions and the 
SMC are much flatter and in the range of 0.63-0.85. 
The slopes tend to follow the thermal fractions of the 
regions which range from 0.5 to 0.95. The thermal frac- 
tion of the RC emission alone can provide the expected 
FIR/RC correlation. The results are consistent with 
a common source for ultraviolet (UV) photons heating 
dust and Cosmic Ray electrons (CRe~s) diffusing away 
from the star forming regions. Since the CRe~s appear 
to escape the SMC so readily, the results here may not 
provide support for coupling between the local gas den- 
sity and the magnetic field intensity. 
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1 Introduction 

One of the most puzzling and robust relationships 
in extragalactic astronomy is the virtually ubiqui- 
tous correlation between far-infrared (FIR) and radio- 
continuum (RC) measurements o f star-forming galaxies 
( Yun. Reddv. and Condon! 200lh (hereafter referred to 



as Y01). The essentially linear r elationshi p extends 



snip 

for over five orders of magnitude (jCondonl Il992t ) ex- 



hibiting the amazingly small scatter of ~ 0.2 dex 
(z dex = 10 z ). Virtually all star forming galaxies 
are included: normal barred and unbarred spiral galax- 
ies, irregular and dwarf galaxies, and Seyferts including 
radio-quiet quasars. This correlation has been found 



to ap ply to galaxies even beyond z=l (jAppleton et al 
20041) . 



The standard explanation for these observations is 
found in models that use the conversion of energy 
resulting from the formation of young massive stars 
(> 20 Mq). FIR emission is posited to be the re- 
sult of dust heated by UV radiation from these mas- 
sive young stars. These stars also provide energy for 
the thermal radio emission found in star forming re- 
gions. Massive stars quickly evolve to the terminus of 
their life cycle and explode into supernovae. The super- 
nova remnants (SNRs) expand rapidly into the local in- 
terstellar medium (ISM) and the shock wave provides 
the mech anism for accelerating Cosmic Ray e lectrons 
(CRe-s) (|Helder et aLll2009tlSchure et adl2009h . Those 
electrons interact with the galactic magnetic field and 
are thus responsible for the non- thermal (synch rotron) 
component of the RC emission ( Condon! 1992 ) . They 
also interact with particles in the ISM and contribute to 
the thermal (free-free) component of the RC emission 
along with UV radiation from hot massive stars. 

This simple picture cannot explain the small disper- 
sion in the FIR/RC ratio measured among the diversity 
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of galaxies which span the observed ranges of star for- 
mation rates (SFRs), magnetic field strengths, metal- 
licity, dust grain chemistry, and ISM mass values. 

The scale at which models break down provides an 
important constraint on the physical mechanisms de- 
tailed in that model. For galaxy-wide scales some mod- 
els work well on spiral galaxies. For very small scales, 
for example, within several kpc of the Sun, the FIR/RC 
ratios of star forming regions are not consistent with 
expected FIR/RC values. Instea d, the radio emission 
appe ars to be primarily thermal ( Haslam and Osborne! 



are accompanied by a res olved region of Population II 
stars with My — —3m. Schwering and Israeli (1993) 



identified 249 infrared sources in the SMC. The mor- 
phology of the Population II stars has been measured 
and appears to suggest a spheroidal distribution with- 
out the prominent irr e gular fe atures seen in the younger 
stars dBolatto et all l2007t ICioni. Habing. and Israel 
2000t iMaragoudaki et all l200ll: IZaritskv et all l2000h 



Earlier studies of the motion of Hi clouds were in - 
terpreted as representing rotation (| Wester fundi 1 1 99 71) . 



1987). 



Hughes et all |2006l ) examined the Large Magellanic 



tion in the SMC ( 


~w — . ~- , — ~ ~ — . - — ~ — 
Piatek. Prvor. and Olszewski 


2008) 


(Costa et al. 


2011). 



Cloud (LMC) for 60 /mi FIR / 21 cm RC correlations 
on spatial scales from 1.5 to 0.05 kpc and found a tight 
correlation on spatial scales above ~ 50 pc. 

The SMC c an be described as a gas-rich late-type 
dwarf galaxy (jBolatto et all 120071 ) . It has a gas-to- 
dust ratio 17 times higher than the Milky Way Galaxy 
( Koornneef and de Boerl 119841 ). It is a member of 
the local group and is classified as irregular (ImlV-V) 
( Sandage et aLlll994l). It may also b e a satellite galaxy 
of the Milky Way dWesterlundlll997l). The SMC is cen- 
tered at 60.6 ±3.8 kpc ( Hilditch. Howarth. and Harries! 
20051) from the Galaxy. The spatial scale of the SMC is 



In this paper, Section [2] discusses the origins and 
characteristics of the data used in terms of the fre- 
quencies of each data set, characteristics of the data, 
and source instruments. Section [3] discusses the anal- 
ysis methods applied to the data sets. Section [4] 
is a discussion of results of the analysis and con- 
siders some FIR/RC galaxy models. Data process- 
ing and analysis for this paper w as done with Miriad 

dSault. Teuben. and Wrightlll995[ ). K arma ( Gooch and Barnes^ 
Il996l). Matl ab dMatlabl l2010l). DS9 djove et all 120031) 
and Ftools ( Blackburn. Payne, and Haves! 1995 ). 



2 Observational Data 



0.3pc/". The line of sight to the SMC has an extinction 
of only AV^0.2 magnitudes and a reddening of E(B- 
V]~0.04 dBolatto et qj2007tlsdilegel. Finkbeiner. and Dav is| 

19981) . The data used for this study are s ummarized in Table CD 



The SMC appears to have been tid ally disrupted by a The F IR data are from the IRAS dMiville-Deschnes and Lagache^ 
recen t close encounter with the LMC ( Murai and Fnjimotnj j 2005h and the Spitzer satellites ( Bolatto et al. 2007 ). 



19801) . T here is no trace of any spiral structure in 



the SMC ^Sandage and Tammannl Il98ll ) but both a 
"bar-like" structu re and a "wing" have been reported 
dWesterlundlll9971 ). The bar structure is probably the 
body of the SMC and contains most of its gas and star 
formation activity. The wing feature continues into the 
bridge which is an H I region exte nding from the SMC 
to the LMC dBolatto et aZ.ll2007h . The SMC may be 
more complex than the LMC considering that various 
population cen troids vary by ~ 2 0' in R.A. and by ~ 1° 
of declination dWesterlundlll997l ). The SMC has an ab- 
solute luminosity of Mb = 2.9m d Westerlundl 1 1 9971 ) . 

The SMC has the lowest m etallicity of any gas - 
rich galaxy in close proximity dBolatto et all 120071 ). 
This makes it one of the best galaxies to provide un- 
derstanding of processes in very high redshift galax- 
ies. Measurements suggest that the dust grains are 
mostly silicates rather than carbonaceous grains which 
distinguish the SMC's ISM fr om that of the Galaxy 
dWeingartner and Draindl2001l ). 

The stars resolved in the SMC are blue and red su- 
pergiants with Mb = —7m. These Population I stars 



The RC data are from the Australian Telescope Com- 
pact Array (ATCA) which consists of five movable 22-m 
antennae, and the Parkes 64m radio telescopes. (For 
reference see the series of paper s on radio-continuum 
studies of t he M ag ellanic Clou d s: lHaynes et all dl 991~): 
" t9921): | Klein et all (fl993 | ): IFilipovic et al 



Xu et al 

— TT 



1995, 1997 




1998 a 



bf) .IFilipovic et al. 



( 2002l);|Payne et alM 



Filipqvic et oil d2005l):lReid et all d2006r).|Pavne et alM 



■ and (jCrawford et alhoiA IWong et aLlboill 



Also, included are data sets for H I dStavelev-Smith et a/.| 
1997: Stanimir ovic et al 1fl999h . Ha dSmith. Leiton. and Pizarrol) .| 



and CO dMizuno et all 120011 ) detailed below. In order 
to make comparisons between the data sets, all of the 
data are regridded to a common pixel size of 5" using a 
gnomonic projection centered at 15° RA and —73° DEC 
(equatorial J2000 coordinates). The whole SMC is de- 
fined as the data shown in Figure [TJ The irregular cov- 
erage of the SMC by the Spitzer data suggested that 
a division of the data sets into several regions would 
be appropriate. The 160 /xm Spitzer data is used as 
an approximate template for the division into the five 
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regions. Figure Q] shows the regional definitions super- 
imposed on all of the data sets as defined by the 160 /itm 
Spitzer data. It is clear that the bar of the SMC is rep- 
resented by regions 1-3 and the wing by regions 4 and 
5. Table [2] lists the coordinates of the boundaries and 
centers of the regions as well as their areas. Table [3] 
contains the flux densities from each region and from 
the whole SMC. 

2.1 Infrared Data 

The IRAS data used in this paper are the im proved re- 



from the center of the ATCA as well as including 
the l ow spatial frequenc y data from the Parkes sur- 



vey ([Havnes et all 119911 ). The addition of this in- 



complete, 2"-3" interference pattern, provides greater 
resolution to the image for the evaluation of previ- 
ously unresolved sources and will be used here to en- 
hance morphological studies using wavelet cross corre- 
lations. The additional data does not cover the full 
extent of the SMC but covers all of the area stud- 



ied in this paper. See lHaynes et all (11991); IXu et al 



dl992l):lKlein et all (119931) : 



1998.a.b), Filipovic et 



proce ssed data from the IRAS satellite survey |Miville-Desc P^QM eM-VKSIOS) 
2005). The data are in two bands: 60 /im and 100 nm 



Filip ovic et 



1995, 



1997, 



20021): ll'avne ( I all I2OO4I): 



Reid et al] ~ (l2006l); IPayne et al 



and have been corrected for zodiacal light, calibration, 
zero levels and striping problems (see Figure [TJ. 

The data used from the Spitzer satellite are from the 
MIPS instrument. They are from the Spitzer Survey of 
the SMC (S3MC) c overing most of the bar and wing 
( Bolatto et all\2007\ ) in two bands, 70 /j,m and 160 nm. 
The data are corrected for zodiacal light, Galactic fore- 
ground and Cos mic Infrared backg round using Leroy's 
published data ( Leroy et all 120071 ) to adjust the final 
calibration. 

2.2 Radio-Continuum Data 



20071) , and lCrawford et all (|201lh : IWong et all (|201lll; 
2012). 



Data from two radio telescopes are used in this study; 
the Parkes Radio Telescope and the Australian Tele- 
scope Compact Array (ATCA) . For high-resolution cor- 
relation investigations, the combined ATCA+Parkes 
data is used and for the low resolutio n thermal frac- 
tion in vestigations, the Parkes data from lFilipovic et all 
( 19971 ) is used exclusively due to its superior overall 
flux-density accuracy. 

The high resolution RC data at 21 cm and 13 cm 
are from observations made with the ATCA telescope 
in 1992 from October 6 th to the 8 th . The observations 
were made using five antennae in the 375 m east-west 
array configuration ( Wong et al. 20121 : Filipovic et < 



20051 ) at 1.42 GHz and at 2.37 GHz. The angular res- 
olution is -98" (FWHM) and -40" (FWHM). There 
were 320 pointings used to create a mosaic of data cov- 
ering an area of (~ 20°) 2 . The "short spacing" low 
spatial frequ ency data fill is accom plished using Parkes 
survey data ( Filipovic et all\1997 ). 

The 6 cm and 3 cm high resolution continuum data 
are from (jPickel et alhoid ). at 4.8 GHz and 8.6 GHz, 
using the EW352 (February 2005) and EW367 (Febru- 
ary and March 2006) configurations in 3564 individ- 
ual pointing positions achieving resolutions of 35" and 
22" respectively. The data were reprocessed to in- 
clude data from a sixth antenna approximately 6 km 



2.3 Hi Data 

The H 1 21 cm data are from observations made with the 
ATCA telescope in 1992 from October 6 th to the U th 
and August 2 nd and 3 rd . The observations were made 
using fiv e antennas in the 375 m east- west array config- 
uration (|Stavelev-Smith et <xZ.|[l995lll997h at 1.42 GHz. 
The angular resolution is ~90" (FWHM). There were 
320 pointings used to create a mosaic of data covering 
an area of (~ 20°) 2 . The "short spacing" low spatial 
frequ ency data fill is accomplis hed using Parkes survey 
data (jStanimirovic et g/J E999h from 1996 March U th 
to 17 th . 

The column densities are correcte d for self-absorpt ionl 



using the correction factors from IStanimirovic et 



( 19991) . The Hi mass is calculated for each region and 
the whole SMC and is shown in Table [3] Results re- 
ported here of 2.8 • 10 8 HI M Q from a 3° x 3° re- 
duced data set compare favorably to a publi s hed v alue 
of 3.8 • 10 8 HI M Q from IStanimirovic et all (|l999h for 
a 5° x 5° image. 

2.4 Ha Data 

The Ha data is from the Magellanic Cloud Emission 
Line Survey (MCELS) p rogram using the UM/CTIO 
Curtis Schmidt telescope (jSmith, Leiton. and Pizarrol) . 
The central 3.5° x 4.5° of the SMC are imaged in 
overlapping 1° x 1° fields with 2048 2 pixels per field 
in such a way as to provide at least two samples of 
each pixel. The Ha filter is centered on 500.7 nm with 
a width of 3 nm. Two continuum channels are also 
imaged, 685 nm, S — 9.5 nm, and 513 nm, 5 = 3 nm 
for stellar contribution subtraction. 

2.5 12 CO j = (1 - 0) Data (115.271 GHz) 

CO data is from the NANTEN sub-millimeter observa- 
tory on Pampa la Bola in the Atacama Desert, Chile 
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Table 1 Data used in this study include IRAS and Spitzer satellite images, Radio images from ATCA and Parkes telescopes, 
Hi images, Ha images from NFCCD/CTIO and CO j = (1 - 0) images. 



Band 


Telescope 


Resolution 


Resolution 


Reference 






(arcsec) 


(PC) 




60 fim 


IRAS 


120 


35.3 


(Mivillc-Dcschncs and Lagachc 2005) 


100 iim 


IRAS 


120 


35.3 


(Miville-Deschnes and Laa'achc 2005) 


70 \im 


Spitzer 


18 


5.3 


(Bolatto et al. 2007) 


160 fim 


Spitzer 


40 


11.8 


(Bolatto et al. 2007) 


3 cm RC 


ATCA+Parkes 


2 


0.6 


(Dickel et al. 2010) a 


6 cm RC 


ATCA+Parkes 


3 


0.9 


(Dickel et al. 2010) a 


6 cm RC 


Parkes 


162 


48.6 


(Filipovic et al. 1997) a 


13 cm RC 


ATCA+Parkes 


60 


17.6 


(Stavelcv-Smith et al. 1997; Stanirnirovic et al. 1999) a 


21 cm RC 


ATCA+Parkes 


90 


26.4 


(Stavelev-Smith et al. 1997; FiliDovic et al. 1997) a 


21 cm RC 


Parkes 


1128 


338 


(Filipovic et al. 1997) a 


Hi 


ATCA+Parkes 


98 


28.8 


(Stavclcv-Smith et al. 1997; Stanirnirovic et al. 1999) a 


Ha 


NFCCD/CTIO 


3 


0.9 


(Smith, Leiton, and Pizarro) 


CO 


NANTEN 4-m 


156 


45.8 


fMizuno et al. 2001) 



Also see lHavnes et aU lll99ll); IXu et all lll992ll; iKlein et all Jl993l); iFilipoyic et all Jl995l. Il997l. Il998llallbh. [Filipovic et all pOO ft : 
Payne et all J2004I' ); iFilipovic et all 1120051 ); iReid et all (120061') ; iPavne et all (120071 1 . and ICrawford et all \201W . IWone et all l l201lllal 



20121) 



Table 2 Coordinates in degrees of the boundaries of the 5 regions as shown in Figure [T] Each region is defined by the 
coordinates of the region's vertices. The centers of each region and the areas are also shown. 



Vertex 


SE 


SW 


NW 


NE 


Center 


Area 


Regions 


RA/DEC 


RA/DEC 


RA/DEC 


RA/DEC 


RA/Dec 


De? 


1 

2 
3 


12.9528/-73.8571 
15.0297/- 73. 1336 
16.5846/-72.5451 


10.4595/-73.5514 
12.4418/-72.8108 
14.0808/-72.2471 


11.7169/-72.7425 
13.3596/- 72. 1624 
14.8865/-71.6050 


14.0390/-73.0425 
15.8403/-72.4741 
17.3130/-71.8711 


12.2920/- 73. 2984 
14.1679/-72.6452 
15.7163/-72.0671 


0.6647 
0.5767 
0.5679 


4 
5 

SMC 


19.4525/-73.7627 
16.8303/-73.3567 
20.5906/-74.4285 


16.6597/-73.5182 
15.2447/-73.1735 
9.4094/-74.4285 


17.6519/-72.6450 
16.0482/-72.4953 
10.2887/-71.4425 


20.2119/-72.8758 
17.6098/-72.6634 
19.7253/-71.4424 


18.4940/-73.2004 
16.4332/-72.9222 
15.0000/-73.0000 


0.7432 
0.3589 
9.0000 



Table 3 Regional measured data and derived results for Hi in terms of M© and H2 mass from CO measurements. The 
percentages refer to the portion of the total summed over all regions for each wavelength shown in this table. The SMC 
data refers to the whole SMC data set which includes data not contained in the regions. 



Data 


Region 1 


Region 2 


Region 3 


Region 4 


Region 5 


SMC 


Units 


60 jim 


24.9(39%) 


13.1(20%) 


16.0(25%) 


7.33(11%) 


2.88(4%) 


92.6 


10 2 


Jy 


100 /im 


47.6(37%) 


26.4(21%) 


29.3(23%) 


17.7(14%) 


7.47(6%) 


208. 


10 2 


Jy 


70 /im 


26.6(34%) 


18.8(24%) 


17.0(22%) 


9.02(12%) 


5.98(8%) 


102. 


10 2 


Jy 


160 fj,m 


73.9(35%) 


42.2(20%) 


45.7(22%) 


34.3(16%) 


15.2(7%) 


332. 


10 2 


Jy 


Hi 


2.64(30%) 


1.87(22%) 


1.41(16%) 


1.79(21%) 


0.954(11%) 


48.2 


10 27 H 1 /cm 2 


Hi 


4.56(30%) 


3.23(22%) 


2.43(16%) 


3.09(21%) 


1.65(11%) 


27.9 


10 7 


Hi M q 


Ha 


44.1(30%) 


30.2(20%) 


54.7(37%) 


15.6(10%) 


4.78(3%) 


149. 


10" 


11 ergs/ cm 2 /s 


CO 


494.(52%) 


121.(13%) 


179.(19%) 


150.(16%) 


4.72(0%) 


82.6 


10 2 


K km/s 


H 2 


71.5(50%) 


13.9(9.7%) 


32.8(23%) 


24.6(17%) 


0.54(0.4%) 


159. 


10 5 


M 


3 cm 


6.28(26%) 


3.80(16%) 


7.61(32%) 


4.49(19%) 


1.57(7%) 


54.8 


Jy 




6 cm 


5.80(29%) 


3.87(19%) 


6.47(32%) 


2.91(15%) 


0.861(4%) 


33.5 


Jy 




13 cm 


7.59(31%) 


5.69(23%) 


7.63(31%) 


2.94(12%) 


0.680(3%) 


35.1 


Jy 




21 cm 


7.55(29%) 


6.46(24%) 


7.59(29%) 


3.69(14%) 


1.17(4%) 


40.0 


Jy 
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20 19 18 17 16 15 14 13 12 11 10 20 19 18 17 16 15 14 13 12 11 10 
RA in Degrees RA in Degrees 



Fig. 1 All of the data sources used i n this paper are shown in this figure. The IRAS FIR data are the improved reprocessed 
data from the IRAS satellite survey (|Mivillc-Dcschncs and Lagachell2005l ). Spitzer FIR data are from the Spitzer Survey of 
the Small Magellanic Cloud (S3MC) ijBolatto et al.ll2007l l. The 160 )im Spitzer data provides the outline for the definition 
of the 5 regions used in this analysis which is shown in the background source panel. The four sources of RC data used in 
this paper are shown in this figure. Hi data from the ATCA and Parkes telescopes are combined to produce the data used 
here. Details of the data are found in Table [l] Data are displayed with North to the top and East to the left and J2000 
coordinates in degrees. These plots are smoothed with a 25" Gaussian kernel and scaled by the square root of the intensity. 
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( Mizuno et all 120011) . The NANTEN telescope is used 
at 2.6 mm wavelength and the observations have a res- 
olution of 2.6'. Regions 1-3 are well covered by this 
data. Region 4 has coverage primarily around N84 and 
there is only incidental coverage of Region 5. 

The conversion factors for the CO data to H2 



mass used in this paper are from iLerov et al\ (|201ll ). 
Leroy calculated 3 different conversion factors for the 
SMC for different parts of the SMC which varied by 
40%. (Leroy's West Region 1; North — !• Region 3; 
East — > Region 4; Average of Leroy's West and North 
— > Region 2 and Region 5). The results of this con- 
version is expected to reveal only the most dense H2 
clouds. The conversion factors applied to the 5 regions 
and the overall average value applied to the entire SMC 
are shown in Table 0J 

2.6 Background sources 

Since the SMC is essentially transparent to 1.42 GHz 
RC radiation and it presents a cross section of several 
deg 2 , many background sources are expected to be vis- 



the SMC RC data. These background sources 


ttalogued by 


Pavne et al. (2004) 


Crawford et al. 


Wone et al. 


(201llla. 20121). This data consist 



of 717 RC sources of which 616 were identified as back- 
ground sources. They are plotted in Figure [1] along with 
the numbers assigned to each of the five regions on an 
image of the 21 cm RC data. 



3 Analysis 

The analysis presented here uses cross correlations be- 
tween several different wavelengths using traditional 
pixel-by-pixel correlations as well as wavelet correla- 
tions. The pixel-by-pixel correlation studies provide 
intensity correlations as a function of position. The 
wavelet cross correlation technique probes morpholog- 
ical correlations which are sensitive to changes over a 
range of spatial scales. FIR/RC correlation is studied 



over the entire SMC as well as the regions defined pre- 
viously. 

All of the correlation analysis with RC data is done 
excluding data at the coordinates of the background 
sources. The mask consisted of Gaussian function with 
a FWHM of 2' applied to the coordinates of each back- 
ground source. The data sets also had a median filter 
of 15" x 15" applied to remove zero pixels and provide 
a slight smoothing. 

3.1 SMC and Regional FIR/RC ratios 

An examination of the FIR/RC correlation in the SMC 
compared to measurements from more distant galaxies 
is shown in Figure [2] This correlation is measured down 
to scale factors of less than 1 kpc by considering each 
of the five regions in the SMC (see Figure [T]). Three 
of these regions cover the bar structure and two cover 
the portion of the wing structure that is closest to the 
bar. The structure of the correlations within the regions 
is shown in images which display the organization of 
the FIR and RC emissions giving rise to the measured 
values. 

Using the 60 [im IRAS IR flux and the 21 cm RC 
flux, a FIR/RC calculation is made for each of the five 
SMC regions and the whole SMC (defined as a 3° x 
3° region centered at RA 15° Dec —73° ). These results 
are directly compared to the survey done by Y01. YOl's 
catalog identified radio counterparts to the IRAS Red- 
shift Survey galaxies with 60 fim IR fluxes of > 2 Jy. 
The catalog includes 1750 galaxies and contains radio 
positions and redshifts. Also tabulated are 1.4 GHz 
radio fluxes from the selected galaxies and their IRAS 
fluxes. 

An ordinary linear least squares fit is performed on 
this data to determine the FIR/RC ratio as defined by 
the data from Y01. The fit equation calculated in this 
paper is Equation [1] 



log(Li.4GHz)- 



\ Hz , 



log 



+ 



' 12.3771V \ 

. Hz ) 



Table 4 CO conversion factors to H2 from ILerov et al\ 
(2011) as assigned to the SMC regions and the calculated 
H2 mass for each region. 



Region 


Conversion Factor 


H 2 




M Q pc- 2 (K km s- 1 )- 1 


*10 5 M Q 


1 


67 


71.5 


2 


53 


13.9 


3 


85 


32.8 


4 


67 


24.6 


5 


53 


0.54 


SMC 


69 


159 



(1) 

For each 5" x 5" pixel of the SMC data sets, the 
value of the FIR/RC ratio is calculated. The sums of 
this ratio for each region is plotted on the graph with 
YOl's data set (Figure [5]). The line described by the fit 
equation is extended through the SMC data sets pre- 
sented in this paper to show how the fitted ratio of the 
Y01 data corresponds with this SMC data. The corre- 
lation is evident in all five regions and for the whole 
SMC. The distance from YOl's regression fit line to 



7 



the regional and complete SMC data are very small: 
from A=0.017 dex to A=0.115 dex. These differences 
changed slightly depending on whether the background 
sources are included or not. The differences for each re- 
gion and the whole SMC are shown in Table [5] including 
the change due to background sources. 

For comparison, Figure [5] also plots the LMC corre- 
lation calcula t ed fro m data published on the LMC by 
Hughes el all (|2006[ ). The distance to YOl's regression 
line for the LMC is calculated here to be A=0.308 dex. 

3.2 Pixel by Pixel analysis 

Pearson's linear correlation coefficient calculation is 
performed between several pairs of data sets, Equa- 
tion [2j For each pair of data sets, the data set with 
the smallest FWHM beamwidth is convolved with a 
Gaussian kernel to smooth it to the same FWHM as 
the lower resolution data. The sampling is then done 
at a spacing corresponding to the beamwidth of the 
lower resolution data to guarantee statistical indepen- 
dence of the samples. In this equation /jv refers to the 
N th image, fjsn refers to the i th pixel of the N th image 
and r p is the Pearson's linear correlation coefficient. 



E(/li-(/l))(/tt-(/2)) 

E(/i*-</i>) 2 E(/ 2 i-</2>) 2 



(2) 



The error in this correlation is calculated as follows: 



(3) 



The formal errors for these data sets are ~10% since 
n, the total number of points in a data set, is ~10 3 when 
the resolution of the data sets are considered. Of course, 
the real error is possibly larger due to calibration errors 
and other systematic errors in the data. 



Table 5 Compariso n of the difference between the data 
fit to data from Yun (|Yun. Reddv. and Condor] 1200 ll ) and 
the current measurements. The contribution to the fit from 
the background sources is shown to contribute less than 10 
percent to the fit discrepancy. 



Regions 


Without BKG 


With BKG 


Difference 




dex 


dex 


% 


1 


0.1159 


0.1153 


-0.52 


2 


0.0431 


0.0433 


0.46 


3 


0.0465 


0.0460 


-1.08 


4 


0.0183 


0.0172 


-6.39 


5 


0.0297 


0.0283 


-4.94 


SMC 


0.0783 


0.0783 


0.00 



Perfectly correlated images would have a coefficient 
of 1 and perfectly anti-correlated images would have a 
coefficient of -1. Table[6]shows various pairs of data sets 
for the 5 regions and the correlations therein. Values 
above 0.7 or 0.8 in Table [5] indicates moderate to strong 
correlation. Those values are printed in bold face type. 

Region 1 and Region 5 have the highest correla- 
tions with the FIR-H I data sets and includes N 19 and 



NGC 267. SMC SNR position data from lFilipovic el al 
(|2005h is used to determine the SNR population of each 
region. Region 1 has the largest concentration of H I for 
the SMC (see Table [3]) and a complex of 5 SNRs around 
SNR B0045-734. There are 4 other SNRs in Region 1. 
Region 5 contains no identified SNRs. Region 2 shows 
very little correlation in any of the data set pairs. There 
are two identified SNRs in Region 2. Region 3 is domi- 
nated by the Hn region N66 and NGC 346, and shows 
good correlations in the FIR-RC data sets as well as the 
Ha correlations with the RC and FIR data. Region 3 
contains 12 SNRs. Region 4 shows moderate to strong 
correlations of the IR-CO data sets for the coolest dust. 
Region 4 is not well covered by the 70 \xm or the CO 
data. There are no known SNRs in Region 4. It does, 
however, contain NGC 456. 

Figure [3] shows the scatter plots of the 21 cm con- 
tinuum and 60 /im IR data for the pixels of all five 
regions and the whole SMC. For each pair of data sets, 
the data set with the smallest FWHM beamwidth is 
convolved with a Gaussian kernel to smooth it to the 
same FWHM as the lower resolution data. The sam- 
pling is then done at a spacing corresponding to the 
beamwidth of the lower resolution data to guarantee 
statistical independence. It also shows two different 
Least Squares Fits to each data set. The blue line is 
the Ordinary Least Squares Fit, OLS and the green line 
is the Weighted Least Squares F it (WLS) calc ulated 
using the Matlab function Iscov ( Matlab 2010l ). The 
OLS assumes the error of each data point is the same 
and consequently weights every data point equally. The 
WLS uses a weighting for each data point which here is 
calculated as the square root of the sum of the squares 
of the values of each axis linearly weighted from 0.01 to 
1 scaled to the value of the data point from the mini- 
mum value to the maximum value of each data set. 

The difference in the slope and offset of these fit lines 
and the averages are shown in Table [7] The last row 
contains the differences in percent of the OLS slope and 
the WLS slope of the fitted lines. It is expected that the 
OLS fit would be a better representation of the entire 
data set whereas the WLS fit is weighted to lower the 
contribution of data with larger calculated errors and 
would describe the strongest data without the effect of 
outlier points. 
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SMC FIR/RC correlation (IRAS and ATCA) With Yun 2001 




5 6 7 8 9 10 11 12 13 14 15 

Log(Luminosity 60um) sols 

Fig. 2 FIR/RC correlations for the whole SMC (defin ed as a 2.5° x 2.5° region center ed at RA 15° Dec -73° ) and the 
regions defined in this paper compared to the data from lYun. Reddv. and Condon! (|200ll ). Regions 1-5 and whole SMC are 
plotted and labeled. (1-plus, 2-circle, 3-asterisk, 4-x ; 5-sq uare, and SMC-diamond). Additionally, the correlation calculated 
from data published on the LMC bv lHuehes et al\ (|2006h is plotted using the star marker. 
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Table 6 Pearson correlations for pairs of data sets. All 
results have a calculated error of ~ 10%. For each pair of 
data sets, the data set with the smallest FWHM beamwidth 
is convolved with a Gaussian kernel to smooth it to the same 
FWHM as the lower resolution data. The sampling is then 
done at a spacing corresponding to the beamwidth of the 
lower resolution data. Values above 0.7 are emphasized as 
having good correlation. 



Regions 



Data Sets 


1 


2 


3 


4 


5 


RC 21 cm 












60 fim 


0.620 


0.686 


0.825 


0.444 


0.263 


70 fim 


0.573 


0.513 


0.585 


0.484 


0.193 


100 fim 


0.574 


0.607 


0.752 


0.410 


0.201 


160 fim 


0.587 


0.458 


0.656 


0.515 


0.168 


RC 13 cm 












60 /im 


0.609 


0.700 


0.815 


0.501 


0.304 


70 fim 


0.610 


0.628 


0.752 


0.720 


0.127 


100 fim 


0.569 


0.638 


0.740 


0.461 


0.244 


160 fim 


0.629 


0.577 


0.661 


0.616 


0.129 


CO 












RC 21 cm 


0.272 


0.059 


0.311 


0.422 


-0.255 


RC 13 cm 


0.296 


0.073 


0.305 


0.347 


0.032 


60 fim 


0.531 


0.130 


0.318 


0.759 


0.116 


70 fim 


0.531 


0.142 


0.405 


0.471 


0.190 


100 fim 


0.554 


0.170 


0.378 


0.769 


0.146 


160 fim 


0.599 


0.168 


0.479 


0.896 


0.181 


Hi 


0.304 


0.205 


0.173 


0.444 


0.293 


Ha 


0.312 


0.018 


0.307 


0.291 


-0.055 


Hi 












RC 21 cm 


0.404 


0.255 


0.059 


0.092 


0.137 


RC 13 cm 


0.443 


0.307 


0.044 


0.137 


0.109 


60 fim 


0.762 


0.584 


0.197 


0.394 


0.748 


70 fim 


0.709 


0.545 


0.209 


0.226 


0.790 


100 fim 


0.810 


0.680 


0.338 


0.472 


0.833 


160 fim 


0.773 


0.708 


0.401 


0.506 


0.857 


Ha 












RC 21 cm 


0.377 


0.659 


0.859 


0.563 


0.329 


RC 13 cm 


0.273 


0.597 


0.830 


0.722 


0.237 


60 fim 


0.205 


0.327 


0.816 


0.516 


0.408 


70 fim 


0.387 


0.290 


0.638 


0.535 


0.018 


100 /im 


0.158 


0.297 


0.731 


0.475 


0.350 


160 fim 


0.438 


0.247 


0.661 


0.179 


0.020 




12 12 

Log(Luminosity 60um) sols 



Fig. 3 Scatter plot of the regions and the whole SMC. 
For each pair of data sets, the data set with the smallest 
FWHM beamwidth is convolved with a Gaussian kernel to 
smooth it to the same FWHM as the lower resolution data. 
The sampling is then done at a spacing corresponding to the 
beamwidth of the lower resolution data. The two lines in 
each plot are the ordinary least squares and weighted least 
squares fits to the data. The OLS is in green, the WLS 
in blue. The contours in the data plot show the density of 
data points and the Gaussian distribution of the background 
noise. Those contours are at 5, 30, 60, and 90% of the peak 
data density. The position of the square is the average of 
the data set. 
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Strong to moderate Pearson correlations in FIR-RC, 
FIR-H I, and IR-Ha in Regions 1 and 3 seen in Table [fj] 
appear to show the effect of two populations of data in 
those regions. The range of OLS slopes in the regions is 
0.594 to 0.978 while the WLS slopes range from 0.638 
to 1.144 over the regions. In every region the OLS slope 
is flatter that the WLS slope. The set of contour lines 
to the left of the figures is from the background noise 
using measurements from the north east corner of the 
whole SMC data set which are not considered to be part 
of the SMC. The measured average value and the RMS 
noise are used to calculate a background Gaussian noise 
distribution and the results of this distribution are used 
to create the noise contours for the background. Also 
shown are four contours which are created by binning 
the luminosity data into a 100 x 100 bin array. The con- 
tours shown are at 5, 30, 60 and 90% of the maximum 
bin count. Visually examining these contours indicates 
that 90% of the data is found with values quite sym- 
metric to the fitted lines, both OLS and WLS. This 
suggests that the WLS is a better representative of the 
data in the regions. The Gaussian noise distribution 
calculated for the measured background noise uses the 
same scaling. It is clear that there is not a significant 
contribution by noise to the scatter plot for any of the 
five regions. However, the whole SMC plot includes 
portions of sky which are not part of the SMC and con- 
tain only noise so that there is considerable overlap seen 
in the lower luminosity values of the scatter plot with 
the noise contours. The plotted square in each data set 
is the location of the average pixel value. 



Table 7 For each pair of data sets, the data set with 
the smallest FWHM beamwidth is convolved with a Gaus- 
sian kernel to smooth it to the same FWHM as the lower 
resolution data. The sampling is then done at a spac- 
ing corresponding to the beamwidth of the lower resolu- 
tion data. Data set averages for RC 21 cm, Log(Luminosity 
1.4 GHz)(W/Hz), and FIR, Log (luminosity 60 fim) sols, and 
fit parameters to y=ax + b using Ordinary Least Squares fit, 
OLS, and Weighted Least Squares fit, WLS. These numeri- 
cal values are shown graphically in Figure[3] The differences 
in slope values are also shown with the WLS always larger 
than the OLS value. For comparison, YOl's data as fitted 
here using OLS gives a = 0.995 and b — 12.38. 



Regions 




1 


2 


3 


4 


5 


SMC 


RC 


13.51 


13.60 


13.50 


13.11 


13.10 


13.07 


FIR 


1.296 


1.169 


1.195 


0.753 


0.762 


0.590 


OLS a 


0.730 


0.594 


0.978 


0.776 


0.725 


0.755 


WLS a 


0.846 


0.638 


1.144 


0.829 


0.752 


0.819 


OLS b 


12.55 


12.86 


12.31 


12.50 


12.56 


12.57 


WLS b 


12.40 


12.82 


12.09 


12.46 


12.55 


12.51 


Aa% 


13.7 


6.8 


14.5 


6.4 


3.6 


7.8 



3.3 SMC FIR/RC ratio map 

In order to relate particular physical structures in the 
SMC with variations in the FIR/RC ratio, Figure [J 
is calculated showing the log ratio of FIR and 21 cm 
RC data. The q value map is calculated using the for- 
mula from IHelou et ail \ 19881) with the IRAS 60 fim 
and 100 [im data sets. The FIR value is calculated in 
equation SJ 



/ FIR 
\Wrn-' 



1.26- 10 



-14 



2-58S , 60^m + <5l00p III 

Jy 



(4) 



The q values of the pixels are calculated with Equa- 
tion [S] 



q = log 



FIR 



3.75 • 10 12 Wm- 2 



log 



Src 



Wm~ 2 Hz- 1 



(5) 



The average value of the ratio is q=2.65. This sug- 
gests that the SMC has a slight FIR excess when com- 
pared to YOl's average figure of q=2.34. Approximately 



Q -73 




+17 +16 +15 +14 +13 +12 +11 
RA in Degrees 



Fig. 4 This figure shows the q value plot of the SMC 
which has been convolved with a Gaussian kernel of 108 pc 
FWHM. Red > 2.95 > yellow > 2.75 > green > 2.55 > 
light blue > 2.35 > dark blue > 2.15 > white. 
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98% of YOl's galaxies are within a factor of 5 of 2.65, 
~1.6 to ~3.0. Figure [4] shows the distribution of q val- 
ues across the SMC from q <2.15 to q >2.95. 

Region 1, with an average q of 2.80, has an a large 
portion of its area with q >2.65. This suggests that 
there is abundant dust in the area. Relatively high 
Pearson correlations of FIR with H I are also seen in Ta- 
ble[S]for Region 1. There are 9 SNRs in this region. The 
strongest RC feat ure in this region i s N 19 containing 



SNR B0045-734 (IPavrie et all 120041) (jCrawford et al. 
201lHWong et aLll201llla l. l2012l ). The other regions have 
extended areas where the RC excess is present. 

Region 2, with an average q of 2.58, has a large por- 
tion of its area with q <2.65. With two SNRs in this 
region, it could be inferred that there is not much star- 
forming activity in this region so that the dust is cooler. 

Region 3, with an average q of 2.70, has 12 SNRs 
with similar areas with q <2.65 and q >2.65. The high 
star formation rate leads to the very high Pearson cor- 
relation seen in all of the IR data maps with the 21 cm 
data. Ha correlates also very well with the RC data 
and the FIR data. 

Region 4, with an average q of 2.75, is dominated 
by background radio source SMC B0109-7330. It has a 
portion of its area with significant dust and gas in the 
eastern part of the region. It is the only region that 
has strong correlation of CO with the IR emission in 
the Pearson correlation results. 

Region 5, with q of 2.73, has very small total energy 
emission, 0.2 - 0.3 of the average of the other 5 regions. 
Scaled for area, it is ^0.4 - ~0.6 of the average. 



The a nalyzing function u sed here is the Pet Hat (PH) 
function (jFrick et aZ.ll200lh . This function is defined in 
terms of its Fourier transform as follows: 



* Ik 



3S 2 (flog 2 ^) :tt< 
= : otherwise 



< 47T 



(7) 



This function specifies an annulus in Fourier space 
with a median radius of 2tt and is non-zero only for the 
specified annulus. 

3.5 Wavelet cross correlation 

Given a two dimensional data set, the energy of the 
transformation as a function of scale factor a over the 
entire plane is: 



M (a) = 



OO />OG 



\W(a,x)\ d£ 



(8) 



oo J — OO 



The wavelet cross correlation coefficient as a function 
of a is calculated as follows: 



r w (a) 



JfWi (a,x)Wj (a,x)dx 

S/Wl («) M 2 (of 



(9) 



where the subscripts refer to the source data sets 
used in the cross correlation. 

The error in the cross correlation is calculated as 
follows: 



3.4 Wavelet analysis 

Wavelets can be used as a tool for scaling analysis 
(jFrick et ali 120011 ). It is a technique that convolves 



the data with a family of self-similar analyzing func- 
tions that depend on scale and location. The family of 
analyzing functions is defined by dilations and transla- 
tions of the analyzing function, also called the mother 
function. In two dimensions, the continuous wavelet 
transformation can be written as follows: 



W(a,x) 



OC (>OC 



— oo J — oo 



X — X 



dx" 



(6) 



In this expression, x = (x, y) and / (x) is the two 
dimensional data set for which the Fourier transform 
exists, (x) is the complex conjugate of the analyzing 
function, a is the scale factor, and k is the normalizing 
parameter. 



Ar w (a) 



(^) 2 -2 



(10) 



where L is the linear size of the data set. For the 
data used here Ar w (a) is < 15%. 

For each pair of data sets used in wavelet cross 
correlation, the data set with the smallest FWHM 
beamwidth is convolved with a Gaussian kernel to 
smooth it to the same FWHM as the lower resolution 
data. The sampling is then done at a spacing corre- 
sponding to the beamwidth of the lower resolution data. 
Figure [5] shows wavelet correlations of several data sets 
with the 21 cm RC in all 5 regions. The smallest scale 
that can be used to analyze the data is the resolution 
of the data. The largest scale is dependent on the over- 
all size of the data set. The smallest linear dimension 
of the data set should be no less than four times the 
largest scale factor used in the analysis. The resolution 
of the 21 cm data is 90" (27 pc) and the finest resolu- 
tion calculated is 80" (24 pc). The correlation is per- 
formed with IRAS 60 jjim and 100 fim data with 120" 
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(36 pc) resolution, Spitzcr 70 [im and 160 fim data with 
18" (5.4 pc) and 40" (12 pc) resolution respectively, Hi 
data with 98" (30 pc) resolution, Ha with 3" (1 pc) 
resolution, and CO data with 156" (46 pc) resolution. 
For the data sets used here, the maximum valid scale 
value for the wavelet cross correlation is 794" (238 pc) 
for Regions 1-4 and 501" (150 pc) for Region 5 except 
for the CO data. 

The CO map is much more irregular than the other 
data sets. Region 1 is close to being fully covered with 
the cross correlation results beginning at 125". Re- 
gion 2 is not well covered and has a large hole in the 
coverage in the center of that region. The only scale 
factors that satisfy validity requirements are at 125" 
and 200", both of which are shown. Region 3 is also 
close to being fully covered with the cross correlation 
results beginning at 125". Region 4 is split into two 
smaller areas. Interpreting the wavelet cross correla- 
tion in the region is problematic. Region 5 is just a 
sliver of data that is too small and irregular to use in a 
wavelet correlation. Wavelet correlations of Regions 4 
and 5 are not calculated. 

3.6 Wavelet cross correlation for FIR with RC 

Table [5] shows the smallest scale factors for which excel- 
lent wavelet cross correlation is measured. Frick et al 



( 20011 ) regarded 0.75 or above as representing an ex- 
cellent correlation between images and will be referred 
to here as the correlation threshold. The scale factors 
are extracted from the data shown in Figure [5] The 
data clearly shows strong morphological correlations of 
dust emission in the FIR with the 21 cm continuous 
emission from the defined regions of the SMC down to 
bcamwidth limited minimum scale factor, ~ 35 pc. In 
Regions 2 and 3, the limit of the correlations is pri- 
marily due to the data resolution. Table [9] shows the 



Table 8 21 cm RC wavelet cross correlation scale factors 
in pc for Region 1 through 5 with data from several other 
data sets. Listed are the minimum scale factors in pc for 
which the wavelet cross correlations are considered excel- 
lent, which is taken to be > 0.75. Note that the BoldFace 
values are limited by the beamwidth of one of the data sets. 



Regions 




1 


2 


3 


4 


5 


IRAS 60 fim 


53.2 


35.3 


35.3 


51.5 


129.8 


IRAS 100 ixm 


57.3 


35.3 


35.3 


55.9 


136.1 


Spitzer 70 fim 


52.4 


28.3 


24.0 


50.6 


50.3 


Spitzer 160 /im 


54.5 


33.4 


24.0 


52.6 


53.2 


Hi 


129.8 




170.1 


88.7 


83.6 


Ha 


30.6 


25.3 


24.0 


45.9 


53.2 


CO 


97.0 




229.0 







o-s 
IRAS60um 










&-& 

IRAS 100 urn 










Spitzer7Dum 










Spitzer160um 








«^ 


HI 






r° 




Halpha 




o^-o o e-e 




y 


CO 











2.0 2.5 3.0 2.0 2.5 3.0 2.0 2.5 3.0 2.0 2.5 3.0 2.0 2.5 3.0 
log(a/arcsec) 



Fig. 5 21 cm RC wavelet cross correlation plots for Region 
1 through 5 with data from several other data sets. The 
data points for Regions 1-4 are calculated at scale factors of 
80", 125", 200", 316", 501", and 794" which correspond to 
24 pc, 37 pc, 60 pc, 95 pc, 150 pc and 238 pc respectively. 
The data set from Region 5 is too small to support wavelet 
cross correlations of 794" so correlations are only shown up 
to 501". Coverage over the map of CO data is not complete 
except for Regions 1 and 3. 
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wavelet correlations extracted from Figure [6] which cor- 
relates the Spitzer IR data with the hig her resolution 
3 cm and 6 cm data, ( Dickel et all 120101 ). The wavelet 
cross correlation with the Spitzer data shows that the 
FIR/RC wavelet correlation scale factors extend down 
to ~15 pc. Within each region, all of the FIR bands 
show similar correlation patterns with RC data with 
respect to scale factor. The coverage of the regions is 
quite good for the Spitzer 160 fim data but less so for 
the Spitzer 70 fj,m data particularly in Regions 1,3, and 
4. There are also some unexpected differences between 
the 3 cm and 6 cm correlations with the FIR data. 

3.7 Wavelet cross correlation for other bands 

HI cross correlation is seen in Figure [5] with 21 cm 
RC in Region 1. The correlation initially dips before 
beginning a slow rise to the correlation threshold at 
^130 pc. Region 2 never achieves correlation threshold 
and Region 3 doesn't until ^170 pc. Regions 4 and 5 
have approximately the same scale factor of ~85 pc. 

Ha correlations from Figure [5] in Table [8] with 21 cm 
RC are essentially at the resolution of the Hi data at 
~25 pc for Regions 1 through 3. Regions 4 and 5 have 
correlation threshold at larger scale factors of ^85 pc. 
High- resolution Ha correlations from Figure [6] in Table 
[9] with 3 cm and 6 cm shows that the Ha wavelet cross 
correlation extends down to much smaller scale factors 
from 9 pc to 33 pc for Regions 1 through 4. 

The CO correlations are very sparse since the map of 
the data is rather incomplete. Region 1 had correlation 
threshold at 97 pc and above. No significant correla- 
tion is detected in Region 2. Region 3 had correlation 
threshold at 229 pc. No other correlations are possible. 



Region 1 Region 2 Region 3 Region 4 



Spitzer70um 










Spitzer"! 60um 










Halpha 











1.5 2.0 3.0 1.5 2.0 3.0 1.5 2.0 3.0 1.5 2.0 3.0 1.5 2.0 3.0 
log(a/arcsec) 



Region 1 Region 2 Region 3 Region A 



Spitzer70um 










pee* 
Spitzer160um 




»</ 






Halpha 











4 Discussion 

The FIR/RC correlation coefficient, q, determined in 
this study is very close to the expected value for nor- 

Table 9 6 cm and 3 cm RC high resolution wavelet cross 
correlation scale factors in pc for region 1 through 5 with 
data from Spitzer and Ha data sets. Listed are the mini- 
mum scale factors in pc for which the wavelet cross correla- 
tions are considered excellent, which is taken to be > 0.75. 



6 cm RC 



Region 


1 


2 


3 


4 


5 


Spitzer 70 /im 


43.2 


19.8 


14.7 


23.1 


137.3 


Spitzer 160 fim 


45.2 


22.7 


17.7 


17.6 


114.1 


Ha 


9.0 


14.0 


10.4 


9.1 


87.5 


3 cm RC 


Spitzer 70 fim 


37.9 


27.6 


16.4 


33.3 


129.4 


Spitzer 160 /im 


39.7 


36.8 


17.6 


35.5 


80.6 


Ha 


12.6 


32.1 


17.10 


33.8 





Fig. 6 6 cm (top) and 3 cm (bottom) wavelet cross correla- 
tion plots for Regions 1 through 5 with the Spitzer FIR data 
sets and Ha data. These graphs extend the cross correlation 
calculations down to 9pc. The data points are calculated at 
scale factors of 31", 40", 50", 63", 80", 125", 200", 316", 
501", and 794" which correspond to 9 pc, 12 pc, 15 pc, 
19 pc, 24 pc, 37 pc, 60 pc, 95 pc, 150 pc and 238 pc respec- 
tively. Region 5 scale factors are limited to a maximum of 
501", 150 pc due to the small size of the region. 
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mal galaxies even though the SMC is structurally dif- 
ferent from the spiral galaxies which constitute most 
of the population on which the correlation is based. 
Spirals generally have a thick underlying non-thermal 
disk. The expected th ermal fraction for the disk of a 
spiral galaxy is ~0.1 ([Condon! Il992i ). The three re- 
gions that cover the SMC bar measured here all have 
high thermal ratios of ^0.54 to ~0.95 (see Figure [7]). 
The existence of the bar, however, is questioned by 
Zaritskv et al. (|2000h . Their argument centers around 
the hydrodynamic interaction between gaseous compo- 
nents providing the physical mechanism for the SMC 
morphology. They further suggest that the SMC is ac- 
tually spheroidal with highly irregular recent star for- 
mation. The overall thermal fraction of the SMC is 
measured here using 1.4 GHz and 4.8 GHz data to have 
a value of ~0.43 which is higher than p reviously pub- 



lished values of 0.15-0.40 for the SMC bv lLoiseau et al. 



(1983) using 1.4 GHz and 2.3 GHz data but lowe r than 
the value of 0.71 measured by Israel et al. ( 2010[) using 
5 GHz and 10 GHz data. However, each of the SMC 
regions measured here except Region 2 has a thermal 
fraction of >0.71. 

The triggers for star formation event s in massive spi- 
rals are believed to be density waves (|Beckll2000l ). In 
the SMC, there is not enough mass to support such a 
mechanism so that the trigger in the SMC is likely to be 
stochastic which leads to a discontinuous formation his- 
tory wit h bursts of star formatio n lasting short periods 
of time (jKTein and Graeve!ll986h . With only two SNRs 
detected in Region 2, there is clearly no evidence of 
a galactic nucleus with its traditional high star forma- 
tion rate. Figure 2] shows the distribution of q values 
in the SMC and there is no indication of the q value 
de creasing from a "nucleu s" to the edge of the galaxy 
as iHelou and Bicavi dl993l) expect in their "leaky box" 
FIR/RC model. Figure [4] does show that most of the 
SMC bar is consistent with q values within YOl's ex- 
pected limits. The dept h of the SMC extends 5-10 kpc 
along the line of sight, (jMartin. Maurice, and Leaueux 
1989). The measured thermal fraction leads to the con- 
clusion that there are not enough CRe~s to generate 
higher non-thermal emissions and/or that the magnetic 
field strength is too low in most regions. The published 
values of the magneti c fi eld are 5—10 u,G along t he bar 



( Ye and Turtle! Il99ll ). ISreekumar and Fichtell (|l99lh 



concluded that the cosmic ray energy density level in 
the SMC is three to five times too low to be in equilib- 
rium with the cosmic r ay sources. This low cosmic ra; 
density is described by iKlein. Grave, and Beckl (|198 
to be a result of the internal pressure from the CRe~s. 
This causes rapid diffusion of the CRe~s out of the 
star-forming regions. This is consistent with the obser- 
vation that the bar has primarily thermal RC emission 



in Region 1 and Region 3 where there is high density 
H I and star formation activity. Region 2 may be popu- 
lated primarily with CRe~s which have diffused out of 
the two star-forming regions in Region 1 and Region 3. 

4.1 FIR/RC Correlation Models 

The FIR/RC correlation models examined here start 
with the assumption that the source of relativistic elec- 
trons and ionizing photons is the supernova rate. UV 
photons from massive stars are absorbed by dust which 
re-emit the energy in the FIR. The relativistic electrons 
from the SNRs interact with magnetic fields or interact 
wi th ioni z ed ga s to produce RC radiation. 



Voelkl ([19891 ) proposed a calorimeter theory which 



assumes that all of the ionizing photons are absorbed 
by dust and all of the energy of the CRe~s are dis- 
sipated by RC emission within the galaxy. This 
does not appear to be consis tent with the SMC since 
Sreekumar and Fichtell (|199ll ) determined that most of 
the CRe~s escape the SMC. 



Ap plying the "leaky box" model of IHelou and Bicav 
(| 1993T) most of the ionizing photons and CRe~s would 
escape the SMC. The observational prediction of the 
model is that q decreases radially from the center of the 
galaxy to the edges of the galaxy. Regions 1 and 3 are 
the main star-forming regions in the SMC at the present 
time. Star formation in Region 2 must be suppressed 
since the Hi content is the lowest of any region except 
for Region 5, see Table [3] Due to the irregular nature 
of the SMC, different position angles from these star- 
forming areas point towards very different environments 
as seen in Figure |31 This does not appear to match 
the gradual decrease in q modeled by Helou and Bicay 
particularly since neither a nucleus nor a disk can be 
id entified in the SMC. 

Hoernes. Berkhuiisen. and Xu ( 19981) (referred to 
hereafter as H01) published a model which decomposes 
both the RC emission and FIR emission into two parts. 
M31 is used as their test case. H01 selection of M31 
relied on four criteria: 

1. M 31 is close to the Galaxy and well resolved in both 
radio and in the FIR. 

2. The radio emission is dominated by the non-thermal 
component. 

3. The FIR is dominated by cool dust emission. 

4. M 31 has a low star formation rate of about an order 
of magnitude below the Galaxy. 

Of these criteria, the SMC is also close enough for 
well resolved FIR and RC measurements. But the SMC 
is not dominated by the non-thermal RC component 
with the overall SMC non-thermal emission at 57% and 
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a minimum non-thermal emission of only 5% in Re 
gion 3, see Figure [7] 

The supernova rate (Rs) for the SMC is about 1/350 
Rssmc = 2.8 ■ 1CT 3 SW • yr' 1 (jFilipovic et al 



years, 



1998) whereas the Galaxy has a supernova r ate of about 
1/50 years, Rs Ga iaxy = 20 • 1Q- 3 SN ■ yr~ x (foiehl et al 



200G). Dividing these rates by the galaxies' respective 
volumes shows that the supernova rate per unit volume 
is higher in the SMC than in the Galaxy: 



Rssmc _ 2.8 ■ 10~ 3 5iV ■ yr- 

VsMC 



15fcpc 3 



Rscaiaxy _ 20 • 10~ 3 SN • yr~ 

VGali 



0.19 



SN 



kpc 3 kyr 



axv 



1200fcpc 3 



0.017 



SN 



kpc 3 kyr 



(11) 



(12) 



Whereas M31 has a star formation rate of an order 
of magnitude below the Galaxy with those two galaxies 
approximately the same volume, the SMC has a volu- 
metric supernova rate an order of magnitude above the 
Galaxy or two orders of magnitude above M31. 

The SMC does not appear to be consistent with the 
criteria needed to apply HOl's decomposition model. 

4.2 Decomposition of the SMC RC into thermal and 
non-thermal emission 

Data from the Parkes single dish telesc ope are used for 



the R C thermal fraction calculations ([Filipovic et 



19971) due to its superior overall flux-density accuracy. 
The Parkes data were converted to a data map with 5" 
resolution using the Miriad regrid program to match 
the 6 cm data set. The thermal fraction of the emis- 
sions at 21 cm are calculated using equation [T3] from 
Niklas. Klein, and Wielebinskil (|1997l ). Se cm and S^icm 
are the measured flux densities of the regions. 



Sftcm r ( ^Gcm 

— Jth 



S2I 



Vile 



+ fr 



V2\c 



(13) 



In equation 1131 Jth is the thermal fraction, f nt h is 
the non-thermal fraction, f nt h = (1 — fth)- The ther- 
mal spectral index used is otth =0.1. The non-thermal 
spectral index is assigned a value of a nt h = 1-09 after 
Bot et al. ( 2010h . The thermal fractions calculated are 
0.81, 0.54, 0.95, 0.84, and 0.76 for Regions 1 through 
5 and 0.43 for the SMC, see Figure [7] The thermal 
fraction calculation assumes that the synchrotron emis- 
sion spectrum can be described by a single power law 
and that the ISM is optically thin. This assumption 
can introduce a bias into the calculated thermal frac- 
tion if the spectral index is not constant. The index is 



expected to steepen in locations that are increasingly 
remote from the area where the CRe~s were initially 
accelerated since the energy loss of CRe~s is ~ E 2 . 
Conversely, in star-forming regions the thermal frac- 
tion will likely be overestimated as the spectral index 
is flatter than in other regions. However, the entire 
SMC as well as each region shows approxima tely 4 to 
9 tim es the expected thermal fraction of 0.1 (jCondon 
1992) which suggests that CRe s are not necessarily 
the source of the RC energy that makes the FIR/RC 
correlation so consistent. 

4.3 Regional Decomposition of the SMC 

The top panel of Figure [7] shows the relationship of 
the slope of the Weighted Least Squares Fit (WLS) of 
21 cm RC verses 60 \im emissions from Table [7] plotted 
against the estimated thermal fractions. The regions 
show a strong tendency towards increasing WLS slope 
with increasing thermal fraction. 

Region 3 has the highest WLS slope, the highest 
thermal fraction, and the highest 21 cm flux density 
of the five regions. It also has the best Pearson corre- 
lation, Table El of any region for RC 21 cm with all 
of the IR data sets, excellent Pearson correlation of 
13 cm RC with the higher temperature IR emissions 
of 60 [im and 70 fxm and excellent Pearson correla- 
tion of Ha with RC 21 cm, 13 cm, IR 60 /jm, 70 fim 
and 160 /im covering cold dust to warm dust emit- 
ters. The locations of the high H 1 concentration cor- 
relates with the IR emissions by way of the high Pear- 
son correlation indicating that the H 1 and dust are co- 
located. It has excellent wavelet cross correlation of 
the 21 cm with 60 fj,m down to the resolution limit 
of 35 pc. Figure [5] shows the wavelet cross correla- 
tion of the higher resolution FIR and RC data (70 fj,m 
and 160 fim Spitzer data with 6 cm and 3 cm RC 
data), indicating that the correlation extends down to 
15pc. This region includes N66 and NGC346 which 
is an Hll region wit h embedded stellar association 
( Bica and Dutra 2000h and the stronge st star- forming 
region in the SMC (jCignoni et al\\mm with 12 SNRs 
identified. 

The next two regions in terms of WLS slope and 
thermal fraction are Regions 1 and 4. They have WLS 
slope values and thermal fractions very similar to each 
other with WLS slope values, ~ 0.85, and thermal frac- 
tions, ~ 0.82. Region 1 contains N 19 and NGC267, 
and Region 4 contains NGC 456 which are H 11 regions 
with embedded stellar associations (jBica and Dutra 
2000). Region 1 has nine SNRs and high Pearson cor- 
relation values for Hi with QOfrni, 100/im, and 160/xto. 
The excellent wavelet cross correlation of the 21 cm 
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Fig. 7 Top Panel: WLS slope of 21 cm RC verses 60 pm 
data plotted against the thermal fraction. The thermal frac- 
tions plotted are 0.81, 0.54, 0.95, 0.84, and 0.76 for Regions 1 
through 5. The SMC estimate is 0.43. Second Panel: Min- 
imum wavelet cross correlation scale factor for 21 cm RC 
with 60 fj,m that had excellent correlation For Regions 1 
through 5, those scale factors are 53.2 pc, 28.3 pc, 24.0 pc 
51.5 pc, and 129.8 pc. The Region 2 and 3 are from the 
higher resolution Spitzer 70 /im data. For the 3° x 3° cen- 
tral portion of the SMC, excellent wavelet correlation is seen 
down to a scale factor of 36 pc. 



RC with 60 fim extends down to 50 pc. Region 4 has 
excellent Pearson correlation of CO with IR emissions 
at 60 fim, 100 fim and 160 pm with the correlation 
increasing as the temperature of the dust emission de- 
creases. Interestingly, the 21 cm flux for Region 1 is 
only slightly below Region 3 but Region 4 has the sec- 
ond lowest flux at about 1/2 of the flux of Region 1. 
The CO data in Region 4 contains NGC 456 but much 
of the remainder of the region is not covered in the CO 
data. 

Region 5 has a WLS slope of 0.72 and a thermal 

fraction of 0.76. The region contain s the open cluster 

NGC 419 which has been examined by Rubele. Kerber. and Gir ardij 
( 2Q10h who suggests that the star formation history of 
the cluster displays continuous star formation over 1.2 
to 1.9 Gyr. This region exhibits no strong Pearson 
correlations in any data set pairs investigated here (Ta- 
ble [6]) . There is excellent wavelet correlation down to 
130 pc which is the largest measured by this study of 
the SMC. Oddly, the Spitzer 70 \xm wavelet correlation 
with 6 cm RC gives a much smaller correlation thresh- 
old, 50 pc, but the Spitzer 70 pm is consistent with the 
60 /im correlation threshold when correlated to 6 cm 
and 3 cm RC data. 

Region 2 has the smallest WLS slope of the regions 
with a value of 0.62. It also displays the smallest ther- 
mal fraction of 0.54. It contains NGC 330 which is an 
open cluster and two SNRs. It contains no strong Pear- 
son correlations. The wavelet correlation threshold is 
at the resolution limit of 35 pc. The higher resolution 
wavelet correlation graphs, Figure [6j show the the cor- 
relation actually goes down to 19 pc. 

The thermal fraction of the SMC is estimated here 
to be 0.43. T his is hi g her t han the typical value of 
~ 0.1 given bv lCondonl (|1992T ) for spiral galaxies. It is 
ve ry close to the value of 0.45 for the LMC measured 
bv lrlughes eToZl (feood) . The average value of the ther- 
mal fraction of the 5 regions is 0.78 over an area of 
2.9 deg 2 whereas the SMC measurement is over 9 deg 2 . 
This suggests that the regions, primarily in the bar, 
emit mostly thermal RC energy but that the larger area 
contains significant non-thermal emission from a large 
population of CRe~s that have diffused out of the bar. 
The star formation activity is in the bar where the H I 
has high density. The non-thermal emission is taking 
place where the Hi has low density outside of the bar. 
The small size of the SMC suggests that the cosmic 
ray leakage into the nearby areas would be substantial. 
The total 21 cm flux from the regions is 26.5 Jy and the 
total flux in the 3° x 3° data set used for these thermal 
fraction calculations is 36 Jy. 

It is interesting to note that the q values that are 
calculated including only the thermal portion of the 
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RC emission satisfies YOl's range of q values for normal 
galaxies for the five regions: 1) 2.90, 2) 2.85, 3) 2.72, 
4) 2.83, 5) 2.85. The "greater" SMC has a q value of 
3.02 which is just at the 98% q value consistent with 
YOl's correlation study. 



5 Conclusions 

We show in this study that the FIR/RC correlations for 
the SMC and for each of the five regions defined here, 
are consistent with YOl's regression from 1750 galaxies. 
The study presented here considers the "greater" SMC 
and the sub-regions of the bar and wing as well. 

Pearson correlation coefficients of RC 21 cm emis- 
sions with 60 /im and 100 /im emissions have positive 
correlation only in Region 3 which is the region with the 
highest star formation rate. Not surprisingly, Region 3 
also has positive Pearson correlation of Ha with RC 
and FIR emissions. Interestingly, Region 1 has twice 
as much Hi and Hi as Region 3 and has the second 
highest star formation rate in the SMC. Its strongest 
Pearson correlations are in terms of H I with FIR emis- 
sions. Since the star formation rate in the SMC is pre- 
dicted to be stochastic, Region 1 may have a suppressed 
supernova rate at this time. 

Pixel by pixel scatter plots indicate two populations 
of correlations in terms of 21 cm RC and 60 /xm emis- 
sions. The highest slope is from Region 3 which also 
has the largest difference between the WLS and OLS 
slope values, 1.14-0.98, 14.5%. Region 1 has the second 
largest slope and the second largest difference of slopes 
of 0.85-0.73, 13.7%. This correlates to the regions with 
the strongest star formation rate having the strongest 
indication of two distinct emission populations. The 
other regions and the SMC have smaller slopes and less 
significant differences of 7.8% to 3.6%. 

The wavelet cross correlation scale factors for the 
FIR with the RC for the regions are all within 20 to 
60 pc except for Region 5 which is «130 pc and may 
have few features at small scale factors. There is no 
clear trend displayed in this correlation. It is interest- 
ing to note that the region with the highest thermal 
fraction, Region 3, and the region with the lowest ther- 
mal fraction, Region 2, have the the smallest wavelet 
correlation scale factors. Ha has wavelet cross correla- 
tion scale factors with the 6 cm RC emission of between 
9 and 14 pc for all of the regions except for Region 5 
which is 87 pc. 

The regions have q values of: 1) 2.80, 2) 2.58, 3) 2.70, 
4) 2.75, 5) 2.73. The measured value of q for the SMC 
is 2.65 which includes portions of the greater SMC that 
are not included in regions comprising the bar and wing. 



The SMC and all of the regions have q values above 
YOl's correlation value of 2.34 but well within the val- 
ues of 1.64 to 3.04 between which 98% of the q values 
for the galaxies measured by Y01 were found. This 
suggests that the RC emission in the SMC is slightly 
deficient compared to the FIR emission. It could be due 
to CRe~s escaping the bar or wing of the SMC with- 
out losing much of their energy. Since the non-thermal 
fraction of the RC emission is so small, the non-thermal 
RC emission is significantly under-luminous considering 
that the non-thermal fraction is expected to be domi- 
nant. This may be due to the small physical size of the 
SMC and a low magnetic field density. 

The SMC displays a thermal fraction that is much 
higher than is common for disk galaxies . Disk galaxies 
typically have a thermal fraction of 0.1 (|Condonlll992l ) 
while the SMC is measured here to have a value of 0.43. 
The regions have even higher measured values of the 
thermal fraction, 1) 0.81, 2) 0.54, 3) 0.95, 4) 0.84, and 
5) 0.76. This suggests that galaxies can have most of 
their RC radiation derived from thermal sources pow- 
ered by the radiation field from high mass stars in H n 
regions rather than by SNRs. There appears to be a 
weak correlation between the thermal fraction and the 
q values. The thermal fraction values above 0.5 corre- 
spond to q values above 2.7 and thermal fraction values 
below about 0.5 have q values of below 2.7. The SMC 
and all five of the defined regions have FIR/RC corre- 
lations consistent with Y01 q values when considering 
only thermal RC radiation. Thermal radiation can be 
the major source of RC emission from a galaxy suggest- 
ing that measurements of the RC emission of galaxies 
could lead to over estimating the super nova rate by 
assuming the 0.1 ratio of thermal to non thermal radi- 
ation. 
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